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Abstract. We present a gauge invariant approach to photoproduction of mesons on nucleons within a chiral 
unitary framework. The interaction kernel for meson-baryon scattering is derived from the chiral effective 
Lagrangian and iterated in a Bethe-Salpeter equation. Within the leading order approximation to the 
interaction kernel, data on kaon photoproduction from SAPHIR, CLAS and CBELS A/TAPS are analyzed 
in the threshold region. The importance of gauge invariance and the precision of various approximations 
in the interaction kernel utilized in earlier works are discussed. 
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1 introduction 

The hadronic spectrum is still the least understood prop- 
erty of QCD. Most theoretical models predict much more 
states than are actually observed so far in experiments, 
see e.g. [1-3]. This is sometimes called the 'missing res- 
onance problem'. The search for missing resonances has 
therefore been an important goal of various experimental 
efforts. The properties of baryon resonances are presently 
under tliorough investigation at several facilities, e.g., at 
ELSA, JLab, MAMI GRAAL, COSY and SPring-8. Due 
to their hadronic decay modes, however, many baryon res- 
onances have large overlapping widths which makes it a 
difficult task to study individual states. In this respect, 
polarization observables can be used as a tool to filter out 
specific resonances in specific reactions. 

A possible explanation of the missing resonances could 
be that these states do not couple strongly to the pion- 
nucleon channels which have provided to a large extent 
the resonance data. A strong coupling of these resonances 
to channels with strange particles could be unraveled in 
photoproduction processes of KA and KS. Such experi- 
ments have been recently undertaken at SAPHIR [4,5] and 
CLAS [6,7] with high precision. More recently, K'^S'^ has 
been measured with the CB/TAPS detector at ELSA [8]. 
These data indicate that resonances so far undetected 
might have been observed, but an unambiguous theoret- 
ical interpretation is still lacking. In addition, the beam 
polarization asymmetry for —>■ K^A, K~^S^ has been 
investigated by the LEPS collaboration at SPring-8 [9]. 

Based on an analysis of cross sections, beam asym- 
metries, and recoil polarizations, the Bonn-Gatchina res- 
onance model, e.g., demands among others the Pii(1840), 



£'i3(1870), £'13(2170) resonances [10,11], whereas in the 
model of [12] a Pi3(1830) is preferred. 

Electroproduction experiments are even more sensitive 
to the structure of the nucleon due to the longitudinal 
coupling of the virtual photon to the nucleon spin and 
might in addition yield some insight into the possible onset 
of perturbative QCD, see e.g. [13,14]. But at low photon 
virtuality experimental data for jp KA, KS are still 
not available. 

Apparently, there are sufficient data to be met by the- 
oretical approaches. In addition to the photoproduction 
data, further experimental constraints are provided by 
pion-induced reactions on the proton. Any theoretical ap- 
proach that aims to describe the photoproduction data, 
ought to be consistent with the corresponding pion-in- 
duced data. 

A successful theoretical approach to meson-baryon scat- 
tering is provided by chiral unitary methods, see e.g. [15- 
20]. In this framework the chiral effective Lagrangian is 
utilized to derive, for example, the interaction kernel in 
a Bethe-Salpeter equation (BSE) which iterates meson- 
baryon rescattering to infinite order. The BSE generates 
resonances dynamically, hence, without their explicit in- 
clusion the importance of resonances can be studied. Chi- 
ral unitary approaches have been implemented quite suc- 
cessfully for photoproduction processes, see e.g. refs. [21- 
24] , but as a simplification only those diagrams were taken 
into account where the photon is absorbed first and then 
the produced meson-baryon pair undergoes final state in- 
teraction. This simplified treatment violates, in general, 
gauge invariance. In order to guarantee gauge invariance, 
diagrams with the photon coupling to any intermediate 
state of the meson-baryon bubble chain must be taken into 
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account [25,26]. One goal of the present work is to study 
the importance of these additional contributions which 
render the amplitude gauge invariant, but have been ne- 
glected in previous works. 

Another simplification in chiral unitary approaches is 
the reduction of the interaction kernel to the on-shell point. 
Although the interaction kernel appears in loops, it has 
been argued that the off-shell components can be absorbed 
by redefining the coupling constants [27]. In this work, we 
do not employ the on-shell approximation but present two 
alternative methods to retain the off-shell components in 
the interaction kernel. For a discussion of these issues, we 
refer the reader to ref. [28]. 

In chiral unitary approaches to meson photoproduc- 
tion, it has been common practice to utilize s-wave pro- 
jections for the meson-baryon scattering kernel and the 
photoproduction multipoles. We will also analyze the ac- 
curacy of this approximation since in our approach some 
higher partial waves are generated through the small com- 
ponents of the Dirac spinors describing the baryon octet 
fields and through the kaon (baryon) pole term in charged 
(neutral) meson production. 

The main goal of the present study is the construction 
of a minimal approach to meson photoproduction based 
on the chiral effective Lagrangian which is exactly uni- 
tary and gauge invariant. The presented method fulfills 
these important requirements from field theory, while at 
the same time any subset of diagrams cannot be omit- 
ted as this would violate unitarity or gauge invariance. 
In this study, we restrict ourselves to the chiral effective 
Lagrangian at leading order. The inclusion of higher chi- 
ral orders in the interaction kernel is straightforward and 
necessary to obtain better agreement with experiment, 
particularly at higher energies away from the respective 
thresholds. This will be the subject of forthcoming work. 

The manuscript is organized as follows. In the next 
section, the effective Lagrangian and the Bethe-Salpeter 
formalism including off-shell components are introduced. 
The gauge invariant extension to photo- and electropro- 
duction processes is discussed in sect.[3l Section|4]contains 
the comparison with experimental data on kaon photopro- 
duction. The phenomenological impact of the violation of 
gauge invariance and the use of the on-shell approximation 
is also discussed. We summarize our findings in sect. [SJ 
while lengthy formulae and a second, alternative method 
for including off-shell pieces in the interaction kernel are 
relegated to the appendices. 

2 Bethe-Salpeter equation 

The chiral effective Lagrangian incorporates symmetries 
and symmetry-breaking patterns of QCD in a model-inde- 
pendent way, in particular chiral symmetry and its explicit 
breaking through the finite quark masses. By expanding 
Green's functions in powers of Goldstone boson masses 
and small momenta a chiral counting scheme can be es- 
tablished. However, the strict perturbative chiral expan- 
sion is only applicable at low energies, and it certainly 
fails in the vicinity of resonances. In this respect, the 



combination of the chiral effective Lagrangian with non- 
perturbative schemes based on coupled channels and the 
Bethe-Salpeter equation (BSE) have proven useful both 
in the purely mesonic and in the meson-baryon sector 
[15-20]. Such approaches extend the range of applicability 
of the chiral effective Lagrangian by implementing exact 
two-body unitarity in a non-perturbative fashion and gen- 
erating resonances dynamically. 

In this work we restrict ourselves to the meson-baryon 
Lagrangian at leading order 

jO% = {mJp,B]-moB)) 

- f (57^75K, B}) - ^{Bj^^,[u^^ B]) , (2.1) 

where the matrix B collects the ground state baryon octet, 
(...) denotes the trace in flavor space, mg is the common 
baryon octet mass in the chiral limit and 

[£1^,5] =a^S + [r'',B] (2.2) 

is the covariant derivative of the baryon field with the 
'chiral connection' 

= 9^ii] - ^ (uty^u + uv'^u'') . (2.3) 

The external vector field is given by = —eQA^ where 
Q = i diag(2, — 1, — 1) is the quark charge matrix. The 
lowest lying octet of pseudoscalar meson fields (j) enters 
the Lagrangian, in matrix form, as 

[/(</)) = u^{(l)) = exp(V2z^) , (2.4) 

JO 

where /o denotes the pseudoscalar decay constant in the 
chiral limit. Furthermore, we use 

V^U ^d^U~i[v^,U], 

and the couphng constants D = 0.8, F = 0.46 [29]. By 
expanding the chiral connection in powers of the meson 
fields, one derives from the effective Lagrangian the lead- 
ing order (fP'BB vertex (the so-called 'Weinberg- To moza- 
wa' (WT) term), which we use as the driving term in our 
Bethe-Salpeter equation. Stated differently, this vertex in- 
sertion is the 'interaction kernel' of the integral equation. 
One finds for the corresponding potential V (which is the 
WT- vertex graph multiplied by i) 

V''^''\i^4,) ^ g''^''\i, + i^) . (2.5) 

Here, qi and q2 are the four-momenta of the incoming 
and the outgoing meson, respectively, and the coupling 
constants are summarized as a matrix in channel space, 
with the entries 

9''^''' = -^{\'^[W\nX'']). (2.6) 
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P = — n- — + ^^p- 

Fig. 1. Graphical iUustration of the BSE for meson-baryon 
scattering. The filled circle represents the full scattering matrix 
and the open square the driving meson-baryon vertex. 



where A° are the generators of the SU{3) Lie- Algebra in 
the physical (particle) basis. In this representation, a dou- 
ble index bj specifies a particular channel consisting of 
a baryon b and a meson j. The baryon and the meson 
propagator, iS and iA, are also summarized as (diagonal) 
matrices in channel space. They are given by 



(2.7) 
(2.8) 



We can now write down the integral equation for the 
meson-baryon scattering amplitude T^i^"-^ in a rather com- 
pact form (suppressing the channel indices, but remember- 
ing the matrix character of the various amplitudes): 

^Vi^^J)zSip-f)Ail)TU,^^;p). (2.9) 

Here, p = pi+qi = P2+Q2 is the overall momentum, where 
Pi and p2 are the four momenta of the incoming and out- 
going baryon, respectively. The BSE is illustrated in fig.[TJ 
Note that we use dimensional regularization throughout. 
The solution of the BSE reads 

T(^2'^i;p) = ^(^2'^i;p) 

+ W{^^,p;p)Gip)[l-W{p,p;p)Gip)]-'W{p,^^:p) 

(2.10) 

with 



9 ,^r ^Mip-m) g. (2.11) 

1 + lug ^ 1 + giu 



W{^^,p;p) can be obtained from eq. (|2.1ip by replacing 
p — m = p, and so on, m is the baryon mass matrix 
with entries 

^bj,at ^ Sbap^^^ ^ ^2.12) 



and the loop integrals are given by 



jbj,ai 

^ ' (27r)' 



■iA''^''''{l). 



(2.13) 



(2n) 



■i[A{l)S{p~ 0]"^^-^ (2.14) 



Note that in eqs. f^UU\i and ^J^, the symbol '1' rep- 
resents the unit matrix in channel space, while matrix- 
valued denominators denote matrix inversion. In order to 
verify that the above expression for T indeed solves the 
BSE, one employs the relation 



(27r)" 



V{<^^J)tS{]/>-t)Ail)W{lf,p) 



= V{^^,p)G{p)W{p, f,p) + W{^^,f,p) - V{^^, i) 

(2.15) 

for arbitrary q. Substituting the expression for T in the 
BSE, eq. (|2.9p and using the last relation, it is straight- 
forward to confirm that T solves the integral equation, as 
claimed. It is important to keep in mind that all the in- 
volved expressions are matrices in channel space, so e.g. 
Im does in general not commute with the coupling matrix 
g. The on-shell substitutions 



'1,2 



1.2 



1,2 



(2.16) 



which occur in the above expressions, provide the connec- 
tion to the alternative method of solving the BSE pre- 
sented in App. [X] The solution given there is completely 
equivalent to the one described above. Here, again, we 
note that p is a matrix in channel space due to the matrix 
character of the baryon mass matrix m. 

Putting the external momenta on their respective mass 
shells, the form of the solution T of eq. (|2.9p simplifies to 



(2.17) 



where we used the shorthand notation W = W{p,p;p) 
and G = G{p). To obtain the complete on-shell meson- 
baryon scattering amplitude, this expression has to be 
sandwiched between baryon spinors u{p2) . ■ ■ u{pi) . Ne- 
glecting the tadpole integrals /m in W , one obtains from 
eq. (I2T7)) 



(2.18) 



which is the form for the scattering amplitude usually 
encountered when utilizing the 'on-shell-approximation' 
(which consists of setting also internal loop momenta in 
the potential V on-shell) to arrive at a simplified version 
of the BSE. The integral equation eq. (|2.9p then reduces 
to an algebraic equation that can be solved by matrix in- 
version to give an expression like above. Retaining 
the off-shell pieces in the exact solution Ton therefore cor- 
responds to keeping the terms proportional to /m in the 
expression for W . 

In [27] it was argued that the on-shell reduction of the 
interaction kernel is justified as the off-shell pieces can be 
absorbed into the vertices. While this procedure simplifies 
the interaction kernel and the treatment of the BSE, it is 
in fact not necessary in order to solve the BSE as explained 
above. We choose not to discard the off-shell pieces and 
evaluate the diagrams in the way dictated by field theory, 
keeping in mind, however, that the off-shell parts are not 
uniquely fixed by any physical requirements and depend 
on the field parametrization [28]. Recall that eq. (|2.4p is 
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not the only possible choice for the parametrization of the 
meson fields, and that the off-shell amplitudes in ChPT 
depend on this choice. 

Another important reason for taking the full off-shell 
dependence of the interaction kernel into account is that 
the implementation of the BSE in a gauge-invariant de- 
scription of photoproduction processes is then straight- 
forward. This is because we evaluate the occurring loop 
graphs according to the rules of field theory which also 
ensures that gauge symmetry can be incorporated in a 
standard way. This is not the case when utilizing the on- 
shell approximation. 

At present, our interaction kernel is restricted to the 
lowest order contact term derived from the effective meson- 
baryon Lagrangian. The extension to more general vertex 
structures in the kernel is possible and will be discussed in 
future work. For now, it is not our aim to achieve perfect 
agreement with experimental data (in which case higher- 
order terms in the kernel would be indispensable), but to 
construct the simplest possible amplitude for kaon pho- 
toproduction that reconciles the framework of chiral uni- 
tary approaches with the fundamental principle of gauge 
invariance in a straightforward way. 

For the sake of completeness, we add some comments 
on unitarity here. From the expression for Ton, eq. (|2.17p . 
one immediately confirms the condition for two-particle 
unitarity 

Im(To-i) = -Im(G). (2.19) 

Our approach to photoproduction presented in the next 
section will be such that unitarity is also satisfied when 
the photon is coupled to the meson-baryon bubble chain 
summed up by means of the BSE, see also [25]. 

The form of the solution Ton guarantees that the above 
unitarity statement holds. However, the careful reader might 
have noticed that the integrals Im and G{p) occurring in 
the amplitude are divergent for c? — > 4 and that the in- 
troduction of appropriate counterterms might spoil the 
simple form of the solution given in eq. p.l7p . To deal 
with this complication one can show by means of a rather 
lengthy calculation that the modification 



in eqs. (|2.10p and (|2.1ip . One observes that Vs has the 
same form as T, with G and Im replaced by —SG and 
—SIm, respectively. The corrections SG and SIm have the 
form of polynomial terms which serve to absorb the di- 
vergences in the loop integrals. The divergences have thus 
been shifted from the loop integrals in T to the new ker- 
nel. Of course, this is not a renormalization scheme in the 
usual sense, since the additional terms 5V in Vs do not 
correspond exactly to counterterms derived from an effec- 
tive Lagrangian which is obvious form the lack of crossing 
symmetry of the amplitude. It has already been noted in 
ref. [28] that the solution of the BSE can not be renormal- 
ized within a usual renormalization scheme. However, the 
foregoing discussion shows at least that altering the loop 
integrals appearing in the solution of the BSE by terms 
SG and SIm is equivalent to certain modifications of the 
potential. In practice, we merely omit the divergences in 
the loop integrals in our expression for T, which will there- 
fore depend on the regularization scale fi showing up in 
the modified loop integrals. 

This form of the off-shell BSE solution T is not partic- 
ularly convenient for the purpose of implementing it into 
the photoproduction amplitude. Therefore, we rewrite it 
by using the following decompositions: 



V ^V + SV = Vs 



(2.20) 



G{p)^Gi{p)p + Go{p), 
W{p,p;p) = Wi{p)p + Wo{p), 



(2.23) 
(2.24) 

where the matrices Go,i follow straightforwardly from the 
explicit expression for the loop integral G which is given 
in eq. (jB.lip of App. [H and 



Wi{p)=g 
Woip) = g 



h 



Mi 



[1 + Imgf 



1 



1 + lyig 



1 



1 



1 + glu ' 
1 



-mg 



^ + lug 1 + glu 



gm. 



(2.25) 



(2.26) 



Note that Gi (p) and Wi (p) are Lorentz scalars which de- 
pend only on the variable p^ . Using these expressions one 
can derive 



of the interaction kernel V in the BSE, where 
Vs{^,,^,;p)^Ws{^^4^;p) ^.^^ 
-Ws{^^,p; p)SG{p) [1 + Ws (p, p-p)SG{p)]-^ Ws (p, ; p) 



l-W{p,p;p)G{p)^Wif + Wo^ 



(2.27) 



(2.21) 



and 



M.g I - gSlu' 



Wi = -iWi{p)Go{p) + Woip)Giip)) , (2.28) 

Wo = 1 - p^Wiip)Gi{p) - Wo{p)Go{p) . (2.29) 
The above results yield 

G(p)[l - Wip,p;p)G{p)]-' = ni{p)f+noip) (2.30) 

1-(5/m5""''^ "^'l-gSlM^" ^ith 

leads to a new integral equation with a solution T,. The _ a,ip)W,-'Wo[p'W, - WoW.-'Wo]-' , (2.31) 

solution T5 IS obtamed from T by replacmg ^y^-/ i ± u i uj , \ , 

G^G-SG " p'Gi(p)[p2VFi - WoW^'Wo]-^ 

_ _ SIm - Go{p)Wr'mp^m - WoW,-'W,]-' . (2.32) 



-SiMif - m)- —g 

1 1 - qSIm 
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All these results can be combined to provide the decom- 
position of the off-shell BSE-solution into the following 
independent Clifford algebra structures: 

+pTjip)+Tsip). (2.33) 
The scalar coefficient functions T„(p) read 

Tab) = LMp)L^ , 

Tab) = [L2^2o{p)+LMp)]Li, 

nip) = Tj{p) , 

nip) = [p'L^niip) + LMp)]Li + , 
n{p) = T^{p) , 

Tj{p) = [p^L^Hiip) + LMp)]L2 

+ [L2n^{p) + L3ni{p)]L^ - glML2 , 
Tsip) = p''[L2f2oip)L2 + LMp)L2 + L2f2^{p)Ll] 



(2.34) 



where the superscript T denotes transposition of channel 
indices and 



Li 



1 



1 + Im9 ' 



L. = -- 



[1 + 5^m]2- 
1 



1 



-gm 



1 + 5/m 1 + /m5 

Note that L\ = L\ and = L2, but in general L\ ^ 
L3. For external on-shell particles, we can make the usual 
substitutions 



'1' ?2 



m = p 



to arrive at the on-shell scattering amplitude 



To, 



on ' on 



(2.35) 



with 



^on' = P^Ti + mTiTO - mT2 - T2m 

-Tarn - mTi + T5 + Te + Tt , (2.36) 
T^°} = /(Ta + T3 + r4 - Tim - mTi) -I- mTam 

-Tsm - mTe -t- Tg . (2.37) 



The functions T„ given in eq. (|2.34p depend only on the 
variable p^ . They will enter the calculation of the vari- 
ous photo- and electroproduction amplitudes in the next 
section. 



3 Photo- and electroproduction 

The Bethe-Salpeter approach discussed in the last sec- 
tion (or the Lippmann-Schwinger equation in the non- 
relativistic framework) can be implemented in electropro- 
duction processes of mesons on the nucleon. In previous 
work, the electromagnetic meson production on the nu- 
cleon was calculated at tree level and the produced meson- 
baryon pair was subject to final-state interactions [16, 24, 
30]. As the photon does not couple to all intermediate 
states of the meson-baryon bubble chain, gauge invariance 
is in general violated and must be restored via artificial 
manipulations. Therefore it seems desirable to develop a 
formalism which implements the principles of gauge invari- 
ance and unitarity in a most natural and straightforward 
manner. We follow here the path which we have already 
outlined in [25], in rather general terms, for the case of a 
photon coupling to a meson-baryon scattering amplitude. 
In this work, we shall be more explicit in evaluating the 
contributions to the various amplitudes in question. Our 
approach for constructing a unitary and gauge invariant 
electroproduction amplitude decomposes into two major 
steps: 

(1) Fix the hadronic part of the amplitude by making use 
of a BSE to implement exact two-body unitarity. 

(2) Couple the photon to the 'hadronic skeleton' construc- 
ted in step (1) wherever possible, i.e. to all external 
and internal lines describing the propagation of the 
involved particles as well as to (momentum-dependent) 
vertices. 

The procedure of step (2) , which is the most natural way 
to guarantee gauge invariance of the amplitude, leads to 
contributions that were usually not considered in chiral 
unitary approaches involving electromagnetic interactions. 
The importance of these additional contributions which 
render the electroproduction amplitude gauge invariant 
can also be quantified within the scheme utilized here. 
The coupling of the photon to internal lines in the bubble 
chain generated by the BSE leads to rather involved ex- 
pressions, since the meson-baryon scattering amplitude T 
appears twice in the corresponding amplitudes, as will be 
made more explicit later when we give the formal expres- 
sions for the various contributions. 

Exact unitarity (in the subspace of meson-baryon states) 
is a fundamental principle satisfied by our approach. It is 
important to note that the procedure of step (2) above 
does not spoil the requirement of unitarity we built in by 
means of the BSE in step (1), as we will show below. First, 
we have to specify the set of graphs which constitute our 
amplitude. 

We start with the tree level B4)B amplitude derived 
from the leading order Lagrangian of eq. (|2.ip . 



(3.1) 



The last index a labels the incoming baryon (here always 
the proton) , while the double index belongs to the outgo- 
ing pair of baryon b and meson i. For fixed a, ^ is a vector 
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Fig. 2. The dressed meson-baryon vertex in the turtle approx- 
imation. The filled circle (open square) denotes the full (tree 
level) meson-baryon interaction. 



in channel space. Furthermore, q is the four-momentum of 
the outgoing meson and 

^ -{X''^{X'\X''})-^^{X''^[X'\X'']). (3.2) 



V2/o 



\/2/o 



To this tree level amplitude, we add the loop contribu- 
tion that accounts for the final-state interaction after the 
meson has left the vertex, to obtain (c/. fig. E]) 



t-p)iS{f~t)A{l)Hg. 



(3.3) 

We call this approximation to the full meson-baryon in- 
teraction the 'turtle approximation', because with some 
imagination the right-most diagram in fig. [2] resembles a 
turtle. Here and in the following we again suppress the 
channel indices for brevity, and p is the momentum of the 
incoming baryon. Using the explicit form of T given in 
eq. (|2.33p . we find 

fj = il^g + T{(^,f- m: p) [{p - m)G{p) - /m]755 

= Wr,{p) + ^r2{p) + pr^ip) + n{p)]j5 (3.4) 

with scalar coefficient functions Fn given by 

Flip) = Tiip^Hi - mHo) + TaliJo - mH^) 

+nHo + TeHi , (3.5) 

r2ip) =g+ P^TiiHo ~mHi) + T^ip'Hi - mHo) 

+p'TiHi + TgHo , (3.6) 

r^ip) = T^iip'Hi - mHo) + niHo - mH^) 

+TrHo + TsHi , (3.7) 

r^p) = p^TsiHo - mHi) + nip'Hi - mHo) 



where 



-p'TrH^ + T^Ho , 



H, = (Goip) - mGi{p))g , 
Ho = {p'Gi{p)-mGo{p)-lM 



(3.8) 



)9- 



The functions Go,i(p) have been defined in eq. (|2.23|) . 

In order to complete step (1), we still have to spec- 
ify which meson-baryon channels contribute in the frame- 
work of our electroproduction model. In this first study, we 
choose to consider only the ground-state octets of mesons 
and baryons, respectively. Moreover, from the topology 
of the hadronic part of the Feynman graph in fig. [2] we 
can conclude that the meson-baryon pairs must have the 
charge and strangeness quantum numbers of the proton. 
This limits the number of channels to six: 



The limitation to these channels can only occur because 
our amplitude is not crossing-symmetric, otherwise more 
channels with different quantum numbers must be consid- 
ered. The violation of crossing-symmetry is a drawback 
of the BSE method which we use to iterate the rescatter- 
ing graphs. To our knowledge, the implementation of both 
crossing symmetry and exact (two-body) unitarity on the 
basis of Feynman diagrams has not yet been accomplished. 
Other methods such as, e.g., an analog of the Roy equa- 
tions (or generalizations thereof) for pion-pion, pion-kaon 
or pion-nucleon scattering, see [31-33], might be required 
to achieve this. Here, our goal is more modest, and we 
sacrifice crossing symmetry in favor of unitarity. 

By now, we have finished the first part (step (1)) of 
our program. Our next task is to couple the photon to 
the hadronic part of the amplitude in a gauge-invariant 
fashion. Inserting the photon coupling at every possible 
place leads to the set of diagrams displayed in fig. [3l 

The first three graphs are generated by coupling the 
photon to the external lines of the hadronic part. This 
leads to the expressions 



Si: = r{^,f)zSif)iier), (3.10) 
SH = {ieQ^Y)^S{f^ - i)r{<^,f^) , (3.11) 
zeQM(2(Z - k)^iA{q - k)r{f - l/^,f^) . (3.12) 



Here and in the following, k is the four-momentum of the 
incoming photon, q is the four-momentum of the outgo- 
ing Goldstone boson and pi and p2 are the four-momenta 
of the incoming and outgoing baryon, respectively, while 
p = pi + k. We have also introduced the meson (M) and 
baryon (B) charge matrices. They are diagonal matrices 
in channel space and read 

Qm = diag(0, 1,0, 1,1,0), 
Qb =diag(l, 0,1, 0,0,1). 

The expressions for the three amplitudes 5*^^^ ^ must be de- 
composed into their independent Lorentz structures which 
is outlined in App. [C] Next we consider the coupling of 
the photon to internal lines of the bubble chain, leading 
to those diagrams that are most tedious to work out. Con- 
sider first the coupling of the photon to an internal baryon 
line, see fig. |3]E. Using the pertinent Feynman rules, the 
contribution of this graph is easily seen to be 



(2^ 



t-p)S{p-t)A{l)eQBr 

xS{p^-t)r{lp^), (3.13) 



pn 



pi], 



AK+, S°K+, S+K' 







(3.9) 



while the coupling of the photon to an internal meson line 
(fig. [3]D) leads to 

J -^Jii^P~Up)S{r)A{p-l)eQM 

X (2(pi -l) + k)'^Aipi - l)r{f^ - t,f^) ■ (3.14) 

Again, the decomposition into independent Lorentz struc- 
tures can be found in App. O 
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F G H 

Fig. 3. Classes of diagrams for kaon production off the nucleon in the turtle approximation utilized here. 



The next class of graphs we consider arises due to the 
'KroU-Ruderman' (KR) term stemming from the covariant 
derivative V^C/ in the chiral Lagrangian. The correspond- 
ing tree level vertex reads (in matrix form) 

'S'KR.tree = T'^TS • (3.15) 

We add the loop contribution that accounts for the final- 
state interaction after the meson has left the KR vertex, 
to arrive at (see fig. [3]F) 

(3.16) 

The remaining graphs are depicted in figs.[3]G,H. These 
contributions arise from the terms proportional to the ex- 
ternal vector field in the chiral connection P^, leading to 
a BcfxjryB vertex rule 

-ier{QM,9}, (3.17) 

and to the following expressions for the Feynman graphs 
of fig. [2]G,H: 

S^Ti^^r{QM,g} J ■^^iS{p^-t)A{i)r{t4^), 

(3.18) 

^WT2 = / ■^J{i,hp)iS{f~l)A{l)Sl^^^ . (3.19) 

The Lorentz structure decomposition of all the above di- 
agrams and the according contributions to the invariant 
amplitudes Bi are given in App. [C] 

Having finished the construction of the electroproduc- 
tion amplitude, we return to the issue of unitarity. The 
crucial observation here is that every electroproduction 
amplitude A^'^(q, k;p) which may be written as (see also 




9Vf ^^{g 9v[g T 

Fig. 4. Illustration of the integral equation (|3.20p satisfied by 
our electroproduction amplitude 



fig. II 

+ / ^,T{q,l-p)iS{f-t)A{l)M'^,{lMp) 

(3.20) 

obeys the requirement of two-body unitarity in the sub- 
space of meson-baryon channels. Here, T is an amplitude 
for meson-baryon scattering that solves a BSE of the type 
of eq. (|2.9p and TMq is a real kernel. 

The proof proceeds in close analogy to the one pre- 
sented in sect. 5 of [25]. Now we note that our electropro- 
duction amplitude given above decomposes into five 'uni- 
tarity classes', each of which obey an integral equation of 
the above type (|3.20[) by construction (with T = T): 

Class 1 : S>^ 
Class 2 ; 5^ + 
Class 3 : 5f -f S"^ 
Class 4 : 

Class 5 : S!^^^ -t- S!^r^2- 

Hence, each class leads to a electroproduction amplitude 
that obeys unitarity for itself, though not gauge invari- 
ance. This is because the expressions ■S'^rpj^, and S!^ 
as well as the tree graphs of S'^j, and 5^ are real in the 
physical region for the electroproduction process and can 
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therefore constitute a kernel A^q in eq. (|3.20p . The sum 

= Si!'+Sli^+S^+S^+S^+S^j^+S^j,j^+S^rj,2 (3-21) 

will then also obey unitarity, due to the linearity of the 
integral equation ()3.20|) in Mq- 

We will now show that the sum Al^ of the five classes 
is gauge invariant by proving k-A4 = for on-shell mesons 
and baryons. This might be obvious from our construction, 
but we include the proof for completeness. The contraction 
of k with the different amplitudes yields 

k- Ss = -er{t^,p) , 

k ■ Su = eQB,r[(^,f^) , 

k-St = eQuri^ - |^,|>^) , 



k- Sb = 



k-S^ 



M 



J-p)iS{f-t)A{l)eC 



M 



/ + $\p)eQMiS{f^ - t) 



(2^) 
X A{l)r{t,f^) 

k ■ Skr = eQugtlb 



k-S^ 



WT2 



{2n) 



^Ti^, l;p)iS{f - I)Ail)eHQM, g} 



dH 



■iS{f^-t)A{i)r{t,f^). 



Here, charged external particles are put on-shell. First, we 
note that the tree graphs are gauge invariant for them- 
selves, since the tree part oi k ■ {Ss + Su + St) equals 

-e-gpb + e-Q^gi^lb + eQugi^ - 1^)75 , 

which is exactly canceled by the first term of fc-SxR (recall 
Qb + Qm — !)■ In order to deal with the loop contribu- 
tions, it is useful to rewrite the integral equations ()2.9|) 
and (|3.3p for T and F, respectively, as 



t]p)^g{i + t) 
d<^l 
{2'kY 



l;p)iS{f-l)A{l)g{t + t). 



(3.22) 



and 



r{t,f,)=gh, + j ^^g{l+tYS{f^-t)A{V)r{t,p^). 

(3.23) 



These equations are equivalent to eqs. p.9p and (|3.3p . Us- 
ing eq. (I3.22P in the second term of fc • S'b and k ■ S'm, 
and eq. (|3.23p in the first term oi k ■ Sb and k ■ Sm, it 
is straightforward to show that the loop contributions in 
k ■ Ai cancel. This, together with the above result for the 
tree contributions, completes the proof of gauge invariance 
for the full amplitude A4. 



4 Results 

In this section, we will perform an overall fit to avail- 
able photoproduction and pion-induced data on the pro- 
ton near the respective thresholds. In more detail, we fit 
the differential cross sections for photoproduction on the 
proton into the K^A, if +17°, K^2J^ final states as well as 
of TT^p K'^A, K^S'^. Inspection of the differential cross 
sections reveals that already at moderate energies away 
from threshold p- and d- waves become increasingly impor- 
tant. Since our approach, which is based on the Weinberg- 
Tomozawa interaction kernel, generates mainly s-waves 
(with one important exception discussed below) and thus 
does not provide a realistic description for higher partial 
waves, we expect it to be valid only in the near-threshold 
regions. We have thus restricted our fits to energy val- 
ues for which the differential cross sections are dominated 
by s-waves, i.e. center-of-mass energies of about 1.80 GeV 
corresponding to photon lab momenta of about 1.25 GeV 
or pion lab momenta of about 1.23 GeV. Still, we will be 
able to extract interesting information from such inves- 
tigations, in particular, we can study in detail the com- 
monly appearing approximations made in the literature as 
already mentioned in the introduction. The extension to 
higher energies requires inclusion of higher order counter 
terms from the effective Lagrangian in the derivation of 
the interaction kernel. However, this is beyond the scope 
of the present work. In this investigation, we focus on the 
importance of gauge invariance and study the on-shell ap- 
proximation in the driving terms and are satisfied with a 
moderate description of both photoproduction and pion- 
induced data on the proton near threshold. Our investi- 
gation sets the stage for systematic improvements in the 
future which will lead to better agreement with experi- 
ment. 

The free parameters in our approach are, on the one 
hand, the three meson decay constants F^, Fk, F,j which 
we vary separately within realistic bounds, as the SU(3) 
symmetry-breaking differences between them are beyond 
our working precision of the effective potential. More pre- 
cisely, SU(3) symmetry breaking is generated by various 
higher order terms in the meson-baryon (meson) Lagrang- 
ian starting at chiral order two (four). Since these contri- 
butions are not included in the leading order WT ker- 
nel, we simulate such effects by allowing variations in the 
various meson decays constants. This, of course, will no 
longer be done when the higher order terms in the inter- 
action kernel have been included. One observes that the 
fitted decay constants tend to larger values reducing the 
strength of the Weinberg- Tomozawa interaction. This is 
consistent with findings in earlier chiral unitary studies. 
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Fig. 5. Differential cross sections for tv p —* A compared to data from [34] (upper panel) and [35] (lower panel). The number 
in each plot denotes the respective cm. energy y^. 



i. e. the WT interaction in many cases produces too strong 
s-waves, see e.g. [24]. We allow Ft^^Fk, and to vary 
between 70MeV and 150 MeV. The values we find for the 
best fits, i.e. the fits with the lowest overall values, are 

F^ = (113... 127) MeV , 
Fk ^ (149... 150) MeV , 

= (74...82)MeV . (4.1) 

Note that Fx tends towards the upper limit of 150 MeV. 
On the other hand, we fit the four different isospin-sym- 
metric scales /i in the loop integrals. For the best fits their 



values are: 

^^jv = (0.46... 0.54) GeV , 
^l^N = (3.29... 4.41) GeV , 
PLKA = (2. 56... 2. 86) GeV , 
/i^^ = (3.66... 4.31) GeV . (4.2) 

We note that these values are roughly in accordance with 
the natural size estimate of ref. [36] (although most of 
them turn out to be somewhat large). 

In figs.iniiniwe present the best fit results for the pion- 
induced differential cross sections tt~p — > K^A, K^E'^, 
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respectively. Although the WT interaction kernel is en- 
tirely s-wave, there are p-vifave contributions due to the 
lower components of the Dirac spinors which are intro- 
duced in the calculation of the scattering matrix T. The 
corresponding total cross sections and the one for tt~p 
K^S^ are shown in fig. [T] We remark that the bands 
in these figures are generated by about 20 fits with a 
very close to its minimal value. Since with the simple WT 
interaction we are not able to describe all these and the 
photon-induced data to a high accuracy, we refrain from 
giving one-sigma error bands as done in our study on K^p 
scattering [37] . This will be done in future work when the 
higher order terms in the interaction kernel will be in- 
cluded. 

The processes 7r~p K^S~ and 7r+p K^S^ are 
dominated by p- and d-waves already close to threshold. 
Hence, a realistic description of these channels cannot be 
provided by the leading WT interaction. In order not to 
overestimate the data we have still included the total cross 
section for 7r~p — s- S~ in our fit. We refrain, however, 
from taking into account data on tt+j? K^S^ at all 
since, regardless of the choice of parameters (within real- 
istic ranges), our approach clearly overshoots the 'K'^p —> 
cross section. This indicates that higher order con- 
tact interactions in the interaction kernel are absolutely 
necessary — particularly for this process — in order to re- 
duce the strength of the WT term. Moreover, it is worth- 
while mentioning that the presented fits exhibit a peak 
in the ■K~p — > rjn cross section due to the S'ii(1535) reso- 
nance. Pion- and photon-induced eta production data will 
be included in future work after the interaction kernel has 
been developed to higher chiral orders. 

Furthermore, we have included in our fit data on dif- 
ferential cross sections of the photoproduction processes 
7P K+A, K+E°, K°Z:+. The resuhs are displayed in 
figs, m m [ini respectively, while the corresponding total 
cross sections are presented in fig. 1111 A few remarks are in 
order: The SAPHIR and CLAS data on charged kaon pho- 
toproduction show some inconsistencies at forward angles. 
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Fig. 7. Total cross sections for tv^ induced production of K'^A 
(top), K^E° (middle), S~ (bottom) as a function of the 
pion lab momentum gi^b- The corresponding cm. energy ^/s is 
also given. The data are taken from [34,35,38] and the compi- 
lation [39]. 



but this can not be resolved within the approximations 
made here. Also, the very different shape of the differen- 
tial cross sections for the K~^A and ivT+Z"" final states can 
be traced back to the isospin selectivity of the A, see also 
ref. [40]. 

To summarize, the results presented here show a rea- 
sonable agreement with data on photon- and pion-induced 
reactions close to threshold, but more realistic interac- 
tion kernels and higher partial waves are required to ob- 
tain better agreement with data, in particular away from 
threshold. This is however beyond the scope of the present 
investigation. 

Of interest here are gauge invariance violations en- 
countered in previous chiral unitary approaches which only 
took a subset of the diagrams in fig. [3] into account. In or- 
der to estimate the typical size of gauge symmetry viola- 
tions, we omit all contributions where the photon couples 
to intermediate and final states and retain only those dia- 
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grams where the photon is initially absorbed (figs.[3]A-|-F). 
Note that both the graphs A and F as well as their sum are 
unitary as explained in detail in sect. [31 For making this 
comparison, we do not refit the parameters but use the 
values obtained in the full approach. Some sample results 
are compared to the full approach in fig. [121 One clearly 
observes that violations of gauge invariance are sizable, 
although this effect could, in principle, be concealed nu- 
merically by readjusting the parameters of the approach. 
This indicates that the additional contributions which ren- 
der the photoproduction amplitudes gauge invariant and 
were omitted in previous work are not negligible and must 
be taken into account. 

In order to be able to compare our results with pre- 
vious chiral unitary approaches we have also worked in 
the approximations employed in these investigations, see 
e.g. [240. To this aim, the interaction kernel for meson- 
baryon scattering is directly sandwiched between Dirac 
spinors and projected onto the s-wave which is then iter- 
ated to infinite order in a geometric series. Furthermore, 
for the photoproduction process only the leading s-wave 



^ Note, however, that in ref. [24] the primary goal was the 
consistent inclusion of the rj' meson in chiral unitary ap- 
proaches. 



(the so-called i?o+ multipole) is considered. Obviously, 
this scheme produces pure s-waves and cannot reproduce 
the structures from higher partial waves in the differential 
cross sections, see fig. 1131 One notes that these approxima- 
tions can indeed be sizable. Of course, in ref. [24] higher 
order terms were included in the i?o+ multipole, but that 
does not change the observations just made. Again, by a 
suitable parameter refitting one might be able to describe 
the total cross sections, but given the more sophisticated 
scheme developed here, such approximations are no longer 
necessary. 

Overall, we have illustrated that both the on-shell ap- 
proximation and the omission of certain classes of dia- 
grams which are required to fulfill gauge invariance con- 
stitute crude approximations utilized within chiral unitary 
approaches in the past. Although these effects can numer- 
ically be concealed to a large extent by readjusting the pa- 
rameters, important theoretical constraints are not met in 
this manner. On the other hand, the approach we have de- 
veloped here is in accordance with these criteria. We stress 
again that our results also show that it is mandatory to 
go beyond the leading WT approximation in the interac- 
tion kernel and to properly include higher partial waves in 
the BSE. Moreover, contributions from three-particle in- 
termediate states, such as irnN, might play a role in kaon 
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photoproduction, see e.g. [41], but are beyond the scope 
of the present investigation and have thus been discarded. 
For a method to incorporate such three-particle cuts, see 
ref. [42]. 



5 Conclusions and outlook 

The search for missing resonances is of great importance 
and actively being pursued at various experimental facil- 
ities. A promising tool to discover new resonances is the 
photoproduction of kaons on protons. A strong coupling 
to channels with open strangeness could reveal new reso- 
nances yet undetected in previous experiments based on 
pion-nucleon physics. Experiments with open strangeness 
are currently performed at ELSA, JLab and at SPring- 
8 providing a host of experimental data. The obtained 
data must be met by theoretical approaches and yield a 
set of tight constraints. In this work, a chiral unitary ap- 
proach based on the combination of the chiral effective 
Lagrangian with a coupled-channels Bethe-Salpeter equa- 
tion is presented. The method is exactly unitary and sat- 
isfies gauge invariance. It improves previous approaches in 
this field which employed only a subset of the Feynman 
diagrams considered here. 

We have fitted both the differential cross sections for 
photoproduction into K^A, K^S''\ K'^S^ as well as of 
the meson-baryon scattering processes Tr~~p K^A, K'^S'^, 
. In the fits, we have restricted ourselves to the 
threshold regions of the respective channels, as we cannot 
expect to obtain a realistic description of these processes 
away from threshold due to the increasing importance of 
higher partial waves absent in the Weinberg- Tomozawa in- 
teraction kernel. In order to improve the agreement with 
experimental data, the inclusion of higher chiral orders in 
the interaction kernel is necessary. In charged kaon photo- 
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production, important contributions to the p and higher 
partial waves are generated by the kaon pole term, which 
is already included at the order we are working. 

Our aim was to estimate, on the one hand, the size 
of gauge invariance violations introduced in previous cou- 
pled-channels analyses by neglecting Feynman diagrams 
where the photon couples to intermediate states of the 
Bethe-Salpeter bubble chain. On the other hand, the sim- 
plification of setting the interaction kernel on-shell and 
performing s-wave projections on the meson-baryon and 
the photon-baryon subsystems is studied. Our investiga- 
tion suggests that both approximations are not justified 
in the treatment of photoproduction processes and lead to 
sizable changes in the results. It is thus important to sat- 
isfy gauge invariance and include off-shell terms in the 
effective potential. The approach presented here fulfills 
these requirements in a minimal way, i.e. any subset of 
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Fig. 12. Sample results which illustrate the importance of contributions which are needed to fulfill gauge invariance. Left panel: 
differential cross sections for photoproduction of K'^A (top), K'^ (bottom) at cm. energies stated in the figures compared to 
data from [7]. Right panel: total cross sections for photoproduction of K^A (top), S'^ (bottom) compared to data from [4]. 
The full result is represented by the solid lines, whereas the dashed lines indicate the contributions from diagrams A+F in fig. |31 
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cos 9 kuh [GeV] 



Fig. 13. Sample results which illustrate the importance of contributions neglected in conventional coupled-channels analyses. 
Left panel: differential cross sections for photoproduction of K'^A at two different cm. energies as stated in the figures. Right 
panel: total cross sections for photoproduction of K'^A (top), K'^ S''^ (bottom). The full result is represented by the solid lines, 
whereas the dashed lines correspond to the outcome of the model [24] given the same values of the fit parameters. The data are 
taken from [4]. 
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the included diagrams cannot be omitted as this would 
violate either exact unitarity or gauge invariance. 

There are three directions in which the approach pre- 
sented here needs to be improved: First, the inclusion of 
higher order terms in the interaction kernel and the re- 
sulting better inclusion of p-waves and higher multipoles 
is straightforward but tedious. This will allow to include 
also the polarization data into the analysis. In this con- 
text, we should also mention that in the chiral unitary 
approach not all resonances are generated dynamically, 
so in fact one might have to include explicit resonance 
fields in certain channels to achieve a precise description. 
Such a method has already been developed for elastic pion- 
nucleon scattering, see ref. [36]. Second, one should also 
include pion and photon induced 77, 77' production data 
as already accomplished in [24] and further constrain the 
scattering and production amplitudes through matching 
to the two-flavor sector (as done e.g. in refs. [43,44]). Last 
but not least, the violation of crossing symmetry due to 
the turtle approximation needs to be repaired. This could, 
in principle, be done by formulating Roy-type equations, 
but is technically involved. 
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A Bethe-Salpeter approach with off-shell 
kernel 

In this appendix, we present an alternative solution of 
the Bethe-Salpeter equation with off-shell pieces in the 
interaction kernel [45,46]. In general, the amplitude for a 
meson-baryon scattering process i?a(Pa)+'/'i(9i) Bb{pb) + 
{qj ) can be written as a function V{t — u, — nic, — 

M|) with masses rric (c = a,b) and Mk {k — for 
the baryons and mesons, respectively, the total momen- 
tum p = Pa + qi = + qj and the Mandelstam variables 
t = {pa—PbY, u = {Pa — qjY- It is uuderstood that the off- 
shell terms 'p^ — uia have been moved to the right, whereas 
the — rrib have been moved to the left. We assume that 
V is analytic in the variable t — u. This is certainly the 
case for all polynomial interactions which are derived from 
contact interactions. The Taylor expansion of y in t — m 
reads 

oc 
1=0 

= Vo + [t ~ u)Vx + [t - uYV2 + . . . , (A.l) 

where we have suppressed the dependence of V and Vi on 
p and the off-shell pieces — nic and pj. — for brevity. 



In the center-of-mass frame the variable t — uis related to 
the scattering angle 9 via 

t-u = 4|pJ|p,| COS0 + 4 {pi - qf) (pI - q]) , (A.2) 

where j| are the moduli of the cm. three- momenta. We 
define a set of orthogonal functions J/ with / = 0, 1, 2, . . . 
(or s,p,d, . . .) which are proportional to the Legendre poly- 
nomials in the center-of-mass scattering angle with the 
first few Ji given by 

Js = l, 

Jp = \Pa\\Pb\ COS^ = - - ^ (pI - if) (pI - I3) ' 

J, = \pj'\Pb\'icOs' - i) = J2 - i |pj2|p^|2 

Jf = |pJ^|Pfc|^(cOs30- ^COS0) ^4-1 \Pa\^\Pb\^Jp- 

(A.3) 

The moduli of the three-momenta in the center-of-mass 
frame can be expressed in terms of Lorentz scalars as 

I |2 1/2,2 2\2 2 

\Pb\'-^{p' + <ij-pir-ql (A.4) 

The expansion of the amplitude V in t ~ u can now be 
reformulated as an expansion in J/, V = J^i^i-^i- o'^" 
der to keep the presentation simple, we restrict ourselves 
to the Weinberg- Tomozawa interaction from now on. The 
generalization to more complicated vertex structures is 
straightforward. For the Weinberg- Tomozawa term V does 
not depend on t — u, 

V = Vo = V, , (A.5) 

and the division of V into on- and off-shell pieces can 
be conveniently accomplished by employing the matrix 
notation (for the reaction ai — > bj) 



^ s 




"'at 1 






rria) , 


^UPb) 


= i^4b- 


mb) , 






2f-ma- mb 
-1 



) , (A.7) 



where g^i^"-"^ is the Weinberg- Tomozawa coupling for the 
channels under consideration as defined in eq. (j2.6p . The 
amplitude V is utilized as the interaction kernel in the 
Bethe-Salpeter equation for the scattering matrix T, c/. 
eq. (l2Jl) . 

+ / ^,V{<^^J)iS{f- t)A{l)T{t4-,p). (A.8) 
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For the decomposition into on- and ofF-shell pieces we 
make the same ansatz for the partial wave decomposition 
of the BSE solution T as for V in eq. ((X6)) . 



rpbj.ai _ T rfibjMi 



(A.9) 



This yields a BSE for Ts which is a matrix equation in the 
combined channel and off-shell space 



(A.IO) 



with a matrix Gs which is diagonal in channel space with 
off-shell submatrices 



G 



ik.ck 



(2^) 



^u,u{t)iS{t)A{p-l)ulS). (A.ll) 



Hence, the integral expression in the BSE factorizes also 
without putting the interaction kernel V and the solution 
T on-shell. In order to solve this equation by matrix inver- 
sion, it is most convenient to decompose the matrices Tg, 
Va , and Gs into a Dirac scalar and a piece proportional to 
(see also sect. [2])- 

B Loop Integrals 



The tadpole integral is given by 



j-bj.ai 

'm 



{27ry P - Af 2 + iQ+ 



n- Mf /M 



gbagjt 



Here is the scale of dimensional regularization and 
A=L(^_ [ln(47r)-fr'(l) + l] 



16772 Vd- 4 2' 

In practice, we will use the AIS renormalization scheme, 
i.e. terms proportional A will be dropped from all ex- 
pressions. In the above result, we have neglected terms 
of 0(4 — d) since we are interested in the limit d ^ 4. The 
diagonal entries in the matrix /m of eq. (|2.13p are given 
by the above expression for the tadpole integral, where 
M is given by the meson mass of the respective chan- 
nel. Moreover, the scale /x can vary between the different 
channels. We shall also define a diagonal matrix /b with 
entries given by tadpole integrals where Mj is replaced by 
the baryon mass rua of each channel. 

In the following, we will leave out the +iO+ prescrip- 
tion in the loop integrals for brevity. The fundamental 
scalar loop integral reads 



loip) = 



(277)-^ [(p-/)2-m2][/2-M2] 



1 



167r2 

M2 



- 1 



4|q| 



2p2 



artanh 



2|q| 



My 



(B.2) 



where 



^(p2 - (to + A/)2)(p2 _ _ M)2) 



(B.3) 



is equal to the modulus of the center-of-mass three mo- 
mentum for a system with particle masses m and M, and 
p2 the total invariant energy squared of the system. 

It is useful to define a diagonal matrix /mb(p), with 
elements given by the above loop integral Iq of eq. (|B.2p . 
where m is the baryon mass ma and M is the meson mass 
Mi of the respective channels. Similarly, we will use ma- 
trices Ibb and /mm, where m and M are both given by 
either the baryon mass rUa (for /bb) or by Mi (for /mm), 
respectively. In the same manner we define 



/ 



bj,ai 
MBB 



{27Tr [{p - I? - m2][(pi - If - mlW - Mf] ' 

An analogous definition applies for /mmb- Only the diag- 
onal elements are needed here, since the coupling of the 
photon does not alter the meson-baryon channel. The loop 
integrals occurring in /mbb and /mmb can be expressed 
in terms of Spence functions [47]. In the course of the 
calculation, loop integrals with a tensor structure in the 
integrand are also encountered. For example, 



d'^l 



{2ttY [{p - If - m^a\\l'^ - Mf 



'mbJ 



(B.5) 



with 



I-'mbJ 



1 

V 



(p' + A/2-m2)/^3"+/^^^«-/^^" 



.(B.6) 

The last expression is simplified for equal masses in the 
propagators, so that, e.g.^ (in matrix notation) 



dH 
(2^ 

We also need 



il^'A{p-l)A{l) 



pi" 



'MM 



(B.7) 



-0^m'^A{p- l)Ail) ^ g^'^ll^lip) 



„2 ^MM\P)-> 



(B. 



where the coefficients of the tensor structures are given by 

{d - l)li!l{p) - (m2 - ^p^)Imm{p) + ^/m , (B.9) 

id - = i^P' - M^)lMMip) + (^ - 1)/m . 

(B.IO) 

Here we have employed the meson mass matrix 

defined in analogy to the baryon mass matrix m, and the 
term p2 the last two equations is of course to be under- 
stood as being multiplied by the identity matrix in channel 
space. 
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The matrix G of eq. (|2.14p can be given in terms of 
scalar loop integrals already defined: 

G{p) = ^ {y~M^ + n?]lM^{p) + lM-I^ +m/MB(p). 

(B.U) 

Integrals with a vector or tensor structure in the numera- 
tor will also be required in the case of three propagators. 
In particular, two of the three masses in the propagators 
will be equal, as already mentioned above. In the follow- 
ing, we will consider the case MBB, noting that the other 
case (MMB) is then given simply by interchanging m and 
M. First, 



{27rr [{p - 1? - ml][{p, - I f - mlW " M|] 

+ S'^'"(pi,p)(pi -p)'' . (B.12) 



The channel matrices A and B defined by this equation 
read 



2D\\ 



1 - 



hlBB + 2/BB(fc) 

"fc2 



-)/mb(pi)- (l + -|-)/MBb)}, (B.13) 



B = ^£<! (4M2 + P- 4p2)/MBB + 2/BB(fc) 



2k^D 



Ap^ 
Ap^ 



/mb(pi) 
/mb(p) 



(B.14) 



Here we have used the following notation: 

f - (B.15) 

M2 ^ - m^), (B.16) 

(B.17) 

k^p-pi, (B.18) 



D = 4p^ - k^ 



fc2 



(B.19) 



and the abbreviation Ap"^ = {Ap'^f. The loop integral 
with three propagators and tensor structure is given by 



jfj.i' 



{2^Y [{p - If - m2][(pi - If - mlW - Mf] 
— 9 

+ C2^^'"(pi +pnp^ +pr + c'j^''\pi-pY{pi-pr 
+ +pY{pi-pr + (pi +pnpi^pT)- 

(B.20) 

The coefficient matrices d read as follows: 
1 



d-2 



1 - Av^ 

MBB + 2^BB(fc) - 2M^A + -^B 



(B.21) 



Co 



k'iM'lMBB+lBBik)) 



+ ^ik'B - Ap'A) - ^(/« (pi) + b)) 

„2 



+ ^(4b(pi) - lilliP)) {d-l)k^C, \ , (B.22) 



4 
1 

2M2((4p2 _ k^)A- Ap^B) 
1 

-4^ 



(4p2 _ k^){M^lMBB + llBBik)) 



-(4p-2- V-A:^)4b(pi) 



(B.23) 



1 



Z\p2 



Ci = ^{ Ap'iM'luBB + lBB{k)) 

-{{Ap^ -e)A- Ap^B) 



+ -(4p-2- V-A:^)4b(pi) 



(B.24) 



A and i? are given in eqs. (|B.13[) and (|B.14[) . while 



MB 



is defined in eq. (|B.6p . Note that in the limit d ^ 4, Ci 
will aquire an additional finite contribution due to the 
divergent terms in the loop integrals. In particular. 



Ci ^ Ci(d = 4) 



1 



(d-l)Ci^3Ci(d = 4) + 



647r2' 
1 



327r2 



(B.25) 
(B.26) 



As already mentioned, the results for the MMB case can 
be obtained from the above formulae by interchanging 
m ^ M. The corresponding coefficients will be denoted 
by A,B and Q. Finally, the preceding results can be used 
to derive the decompositions 



d'^l 
(2nY 



t)eQBrS{f^ - t)A{l) 



+ pf^Fe + pt'pFr + p^j^iFs + p^^fFg + ptn^iF^ 



and 



dH 
" / (2^ 



H = 



(B.27) 



iA{p - l)eQM{p + Pi- 2irA{p^ - l)SU) 



= l''Fi 



1^F2 + i^j^iFi +fi^'pliFi+p^'F^ 
+ p'tFe+ p^'ppT + p^j^iFs + ptfFg + pt^iFw , 



(B.28) 
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where the coefficients of the Lorentz structures are given with 

by 



e 



Ss — — ^ — (r2(p) ~ m„ri(p)) , 

+ pI{A + B)+p'{A-B)), sV = - m^r^ip)) 

F2 = cQbtuImbb , 



mi — s 



t2 

P 

e 



F4 = egB(/MBB - 2 A) , 



Fe = -2eQBm{B + A) , 

F7 = 2eQB{B - A + C2 + C3 - 2C4) , Here, nip is the proton mass and s = p^ ^ {pi + kf . 

V o n (n n \ Moreover, 

^9 = ^8 , = q^^^^Sl + ^7^75^^ + ri^Sl (C.2) 

i^io = 2eQB(C2 + C3 + 2C^-A-B), 

Fi = -2eQMC'i , 

F2 = F3 = F4 = 0, 52 = -25r = ^^r3, 

F5 = eQM(m/MMB + 2m(B-i)), ^Qb . . 2^ ^ , 2^ 

/r^Jnn SZ = —^{ri{p^)mp{u-mj,)-r2{pi)(u-m^) 
Fq = eQM[mlMMB - 2m[B + A)) , m - u 

F 



with 



7 



eQuiA ~B- 2{C2 + C3 - 2C4)) , + K " "^p"^)^3(Pi) + (m - , 

Fg = egM(i + S-2(C'2-C'3)), and 

Fg = egM(i - B - 2(C2 - C3)) , = (p^^f + rfSf^ + + ^q^sf + fq'^sf' 

Fio = eQM (i + B - 2(C2 + C3 + 2C4)) . _^ ^pM^^^P ^ ^pP^^^P _^ ^pM^^^Pi + ^p^^gtP^ ^ 

The following relations are helpful, e.g. when checking (C.3) 
gauge invariance: 

Gi(p) = p2^4 - ^(i^4 + F7 + Fg) 



with 



,2 





5f = 


t - AP ' ' 




-sf, St 






2eQM ,^ 
i - M2 ' ' 






A, 






givi _ 





'2Sl, 

■{Fq-F-;)-Fi 2eQM ,^ afi oil oiP 



e\ . Ap^ - -'T^jp^^' ^* ' -J^NP^^' 



Gi(pi) p^F4 H (F4 + Fg + Fio) The foregoing results deserve some comments. We use the 



k 



2 



abbreviations 



2^^' ri=™2^i(pi)~mp(F2(pi) + F3(pi))+F4(pi), (C.4) 



Ys = (rrip + mpTO)Fi(pi) - (m + TOp)F2(pi) 



and also 

fc2 Zip2 



-"ipF3(pi)+F4(pi). (C.6) 

Go(p) - Go(pi) = y (-^6 - ^5) ^(2^2 + F5 + Fg) Furthermore, we have anticipated that the full amplitude 
;^2 ^ 2 will be multiplied by spinors, so we set 1*^^75 — > —rrip 75. 

The Mandelstam variables u and i are understood to be 
diagonal matrices in channel space and are defined via the 
center-of-mass (cm.) relations 



2 ' " 2 

As a further check of the calculation we have utilized the 
routines provided by the FeynCalc package [48] . 



t = aP + fc2 - 2FfcFg + 2|k||q| cos0, (C.7) 

C Decomposition of the amplitudes u + ml + ni^ + Ap - s - 1 , (C.8) 

We start with the amplitudes of eqs. ((XTU1) - ([XT^ . where |q| and |k| are the moduh of the cm. three-momenta 

of the outgoing meson and the incoming photon, respec- 
Si^ = {(^fSf^ + (^SV + fSV + S'2')7''75 (C.l) tively (see eq. dHSI), B is the scattering angle in the cm. 
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frame, and the cm. energies are given by 



Eg = V|q|2 + M2. 



Of course, multiplication of scalar quantities (such as fc^ or 
rUp) by the unit matrix in channel space is implied where 
necessary and fractions involving matrix-valued denomi- 
nators denote matrix inversions. 

Now we come to the graphs where the photon couples 
to internal meson or baryon lines. The decomposition of 
reads 



(C.9) 

with 



c7 



TUplTsFs - T^mFa - T^eQBGiipi) 
+ sinFa + T^Fi + TrFi - T3mF4)]Y3 

- TgFiFs + TsmFiFs + TscQbGoCpOFs 

- neQBhiY2 + TsGa{p)eQBY2 - nmGo{p)eQBY2 

- siTsFiYs + T5F2Y3 + TJF2Y3 - T3Goip)eQBY2 

- T3mF2Y3 ~ T5Gi{p)eQBY2 

- TjGi{p)eQBY2 + T3mGi{p)eQBY2) , 
= mpinFsYa - TsmFg^s + sUFsYa] 

- TSF5Y3 + T^mFsYs 

- siT3F5Y3 + T^FrY3 + TjFjYs ~ TgmTO) , 
mp[TsFwY3 - T5mFwY3 + ST3FWY3] - Tgi^gia 

+ TsmFeFs - sCTsFg + T5F9 + T7F9 - T3mF9)Y3 , 

5^'' = mp[TeF3 - TamFs - TsegeGiCpi) 
+ s(TiF3 + T2F4 + T4F4 - TimF4)]Y3 

- TeFiYa + TamFiFs + T2eQBGo{pi)Y3 

- T2eQBlMY2 + TeGoip)eQBY2 - T2mGoip)eQBY2 

- s{TiFiY3 + T2F2Y3 + T^F2Y3 - TiGo{p)eQBY2 

- TimFaFs - T2Gi(p)eQBl2 - T4Gi(p)eQBl2 
+ TimGi{p)eQBY2), 



S'b 



- mp[T5F3 + TrF3 + T^F^ - T^mFs 

- TseQBGiipi) - T^mFi + sT3F^\Y3 

- T5F1Y3 - Tri^iFa - ^8^^2115 + Tami^iFs 
+ T3eQBGo{pi)Y3 - T3eQBluY2 

+ T5Go{p)eQBY2 + TsmFaFg + TrGo{p)eQBY2 
+ TsGiip)eQBY2 - T3mGo{p)eQBY2 

- nmGi{p)eQBY2 + s(r3Gi(p)eQB>2 - 13^^213) , 



Si^^ = mp[T2F3 + TiFs + T^F^ - TimFg 

- TieQeGibi) - T2m.Fi + .sTi^ilFa 

- T2F1Y3 - TiFiY3 - Te,F2Y3 

+ TimFiY3 + TieQBGo{pi)Y3 - TieQBlyiY2 
+ T2Go{p)eQBY2 + TamJ^sFa + TiGo{p)eQBY2 
+ TeGi{p)eQBY2 - T^mGo{p)eQBY2 

- T2mGi{p)eQBY2 - s{TtF2Y3 - T^G^{p)eQBY2) , 



,[T(,Ff,-T2m.Ff, + sT^Fs]Y3 



- TeFgFg + T2mFr,Y3 

- s{TiF5 + T2F^ + T^F^ - TimFj)Y3 , 

5f = mp[T5Fs + TrFs - TamFslFa 

- (T5F5 + T7F5 + TsFr - TsmFg 

- TgrnFy + sT3Fr)Y3 , 

Si^" = mp[T2Fs + T^Fs - T,mFs]Y3 

- (T2F5 + TiF^ + TeFr - T^mF^ 
-T2mF-j + sTiFj)Y3, 

sl"' = mp[TeFw - T2mFw + sTiFw]Y3 

- TqFgYs + T2mFeY3 

- s{TiF6 + T2F9 + T4F9 - TimFQ)Y3 , 

S^' = mplnFw + T7F10 - TsmFiojys 

- (T5F6 + T7F6 + TsFg - TsmFe 

- TsmFg + sT3Fq)Y3 , 

Si^"' = mp[T2Fw + TiFw - TimFw]Y3 

- TamFg + sT^F^)Y3 , 



while the decomposition of reads 



s'iA = b^s^^ +p'^si, +pts^^ + ^rst] +fr''S] 



M 



with 



= nmF,Y3 - TsFiYs - neQMdiY2 - ST3F1Y3 , 
5^ = mp[TsFsY3 - Tsmi^g^s - T^eQMd2Y2 + ST3FSY3] 

- TSF5Y3 + nmF5Y3 - siTsF^Ys + nFrYs 

+ TrFrY3 - T3mFrY3 + T3eQMd2Y2) , 
Sm = mp[TsFioY3 - TsmFioYs + neQMd2Y2 

+ ST3F10F3] - 78^6^3 + TsmFeFs 

^ s{T3Fq + T^Fg + TrPg - T3mFg)Y3 

+ sT3eQMd2Y2 , 



T2mFiY3 - TeFiYs - sT^FiYs - T2eQMdiY2 , 



= T3mFiY3 - T^FiYs - TjFiYs - T3eQudiY2 , 



Timi^iFs - T2F1Y3 - T^FtY3 - TieQMdiY2 , 
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M 



M 



- T2eQMd2Y2] - T6F5Y3 + TsmFgla 

- s(TiF5 + T2F7 + T^Ft - TimFj)Y3 

- sTieQMd2Y2 , 

rup [T5FSY3 + TrFsY3 - TgmFsFa - T3eQMd2>1i] 
(T5F5 + ^7^5 + W - TsmFs - T5mF7)Y3 
sTsFrY^ - T5eQMd2Y2 , 

mp[T2FsY3 + T^FsYs - TnnFsYs - TieQMd2Y2] 
{T2F5 + TiF5 + TQF^ - TimFs - T2mF^) 
sT^FrYs - T2eQMd2Y2, 
mp[TeFioY3 - T2mFioY3 + sTiFwYs 
T2eQMd2Y2] - TqFqY^ + Tzmi^g^s 
sTieQMd2Y2 

s{TiFe + T2F9 + T4F9 - TimFQ)Y3 , 
mplT^FioY^ + TrPioYs 

(T5#6 + Type + nPg - T^mPe - 
sT3p9Y3 + T5eQMd2Y2, 
mp[T2pioY3 + T4pioY3 

- TimFioYs + TieQMd2Y2] 

- {T2pe + TipQ + npQ - TimPe ^ 

- sTiFsFg + T2eQMd2Y2 . 
Here we have used the abbreviations T!, 



'-'1 



M 



T5mFg)Y. 



S 



M 



T2mFg)Y. 



di = -2 
1 



- ifc2)jMM(fc) + ^/m 



T,{p) and 
+ S , 



d2 
5 



luuik) 



1 



2 

- 



(fc^-Af^)/MM(fc)+/l 



M 



The Ti{p) are defined in eq. (|2.33p . while the Fi and Pi are 
defined in eq. (jB.27p and (jB.28p . respectively. The term 5 
stems from the limit d — > 4 in eqs. (|B.9p and (jB.lOp . 

Next we turn to the class of diagrams derived from the 
'KroU-Ruderman'-term (called 'Class 4' in sect. [3]): 

^^^R = (s'kk + A +A + iptllh'!^ (c.ii) 

with 

Skr. = eQug + [{Ts - T^m + sT'ijGoip) 

+ s{Tr-, + Tt - T3m)Gi{p) - T^IujeQug , 

41 = 1(^6- Tarn + sri)Go(p) 

+ s{T2 + T4 - Tim)Gi{p) - T2lM]eQM9 , 

st^R = m+Tj-nm)G„{p) 

+ (sTs + Ts - Tz'm)Gi{p) - T'ilM\eQyig , 



stll = [{T2+T,-T,m)Go{p) 

+ (sTi + Te - T2m)Gi(p) - TiIu]eQMg ■ 

Here, the first term in the first line is the contribution from 
the tree graph, and Ti = Ti{p). The last class of graphs is 
^WTi + ^WT2' where 



WTl 



(C.12) 



with 

^WTi = {Qyig + 9Qm){{Gq{pi) - mpGi(jpi))Y3 - IMY2) , 
and 

SwT2 = ('5'wT2 + ^'S'v?T2+^'S'w't2 + ^^'5'^T2)7^75 (C.13) 

with 

SlvT2 = HT^Goip) + (T5 + T7 - T3m)Gi(p)) 
+ {Ts - T5m)Go{p) - T5Im]S'^j.^ , 

sVt2 = HTiG„{p) + {T2 + T4 - Tim)Gi(p)) 
+ (7(3 - T2m)G^{p) - T2Im]Swti > 

^WT2 = [(^^3 - T^m + Ts)Gi{p) 

+ {T5 +Tj- T3ni)Ga{p) - T^ImiS^j,^ , 

st^l^^[{sTi~T2m + n)Gi{p) 

+ {T2 +T4~ Tim)Go{p) - TiIm\S^j,^ . 



Adding the contributions of eqs. (jC.ip - (|C.3P and eqs. (jC.9P - 
(jC.lSp . we can decompose the full photoproduction am- 
plitude (see eq. (|3.21l) ') as follows: 

= 7^75X1 + (7^75Al2 +^^75X3 +^^75X4 

+ ^7^75X5 + ^7''75A^6 + (^fl^l5M7 
VJg -H^p^7,57Wio 
^^lbMi2+ptlbMi3 



+ fp^l5Mii - 



(C.14) 

This can be simplified, using the Dirac equation and mo- 
mentum conservation, to arrive at the operator basis given 
by the Afj^ co mmonly used for photoproduction processes, 
see eq. (jP.ip . The relation between the corresponding co- 
efficients Bk and the functions Aij used in eq. (IC.14P is 
given in eq. (|D.2p . The decomposition of the amplitudes 
into the various Dirac structures is obtained by employing 
the FeynCalc package [48]. 



D Invariant amplitudes 

We consider the reaction 

j{k) + p{pi,mi) B(p2,m2) + M{q,M^) 

and define the Mandelstam variables as usual, 

s = {pi + kf , 
u = (pi - qf , 
t^{p2-pif. 
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They obey the constraint 

S + t + U = m\ + 777,2 + ^'^l + ■ 

The amphtude can be decomposed as 

8 



(D.l) 



where the operator basis is given by 



Here, P — \{pi +^2)- The relation to the coefficient func- 
tions M.j used in eq. (jC.14[) is 



B2 = ^Als + ^-^4 + M5 + Me - m-iMj 

- + m-2)A^10 - ^miMll + i(s + 777 l 7772)All2 

- + m2)Mi3 - ^niiMiA + ^(s + 77717772 )A1 15 , 
B3 = + + \Mi + + 

- ^7772A^7 " + m2)M8 - ^miMg 

- ^{mi + m2)Mia - ^miMn + ^(s + 77717772 12 

- ^{mi + m2)Mi3 - ^miMiA + ^(s + 77717772 15 , 
Bi = ^Alg - ^Mi + + ^Me - ^m2M7 

- -^{mi + m2)Mio - -niiMn + -(s + 77717772 12 

-I- ^{mi + m2)Mi3 + ^miMu - ^(s -I- 77717772)7^15 , 
B5 = -A^i + (7772 - mi)M5 - miMfi 

- (s - 777 l 7772)A^7 , 

Be = -A/Jio - Mil + r772A^i2 

- AI13 - Mu + r7727Wi5 , 

B7 = -^AllO - ^-^11 + ^™2Xl2 

+ ^>^13 + ^-^14 - ^rn2Mi5 , 
Bs = -Ms -M9- ^Mio - ^Xii + ^■m2Mi2 



^Mi3 - ^Mi4 + ^m2Mi5 ■ 



(D.2) 



For a gauge invariant amphtude, the following rela- 
tions for the Bi hold: 

k^Bi + 2{k ■ P)B2 + 2{k ■ q)B3 + 2fc2B4 = 0, 
B5 + {k- P)Be + k^Bj + [k ■ q)Bs - 0. 



Given these relations, one can eliminate two of the Bi 
(conventionally, one takes B3 and B^) and rewrite the 
amplitude in a manifestly gauge invariant form: 



Tfi = iu2 A.MiUi 

i=l 



(D.3) 



with the operator structures 



1 



Ml = ^151^.1.F^^'' , 
M2 = 275P4(Z - ifc).i^''^ 
M3 = 757m9--F^''" , 

M4 = 2757^P,F'^"^ - (7771 + r772)Mi , 

Me = 75fcM7.J^''" , 

where P'^" — e'^k" — e^k^. The particular form of M4 has 
been chosen such that 



iu2M4Ui 2yjm\m2 q • (k x e) 

in the nonrelativistic limit, where the baryon masses are 
large compared to the meson masses and three-momenta. 
The Ai are related to the Bi via 

Ai^ Bi - 2 (™i + ™2)P6 , 



A 



Ml - 1 

A3 = ~Bs , 
1 



B2 



A4 
^5 



7,7 — 777T — 777S 



^\Bi- 
Ae^By. 



s — u + m% 



-Bo -I- 2Bi 



Following the conventions of Chew, Goldberger, Low and 
Nanibu (CGLN) [49], and Berends, Donnachie and Weaver 
[50] , we rewrite this once more, making use of the standard 
representation of spinors and Gamma matrices. In terms 
of Pauli spinors and matrices, one obtains 



1 



6 

^ l — l 



(D.4) 



where the matrix F reads 

F = 7(7 • b J^i + (T • q (T • (k X b) J^2 + icr • k q ■ b JF3 
+ 7(7 • q q • b — 7(T • q 60-^7 — icr • k b^^Tg, . 



Here, the hat over the three-vectors of course means nor- 
malizatit 
fined as 



malization to a unit vector, and the four-vector is de- 



b -e -—k 
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By substituting the standard representation of the Dirac 
spinors and matrices on the left-hand side of eq. (jD.4|) , one 
finds the following expressions for the so-called CGLN- 
amplitudes J-f. 



Ax 



k ■ q 



{yG ~ni2 y^—^ )A4 



^ -mi 



A. 



's — mi 



-^2 = (Vs + ITT-i) 

-Al 



^ -mi 

q||k| 



A, 



S-KyG {El +mi){E2 + m2) 
k ■ q 



y/s + mi 



A. 



, r- k ■ q . 
+ {^/s + m2 )A4 - 



fc2 



= {Vs + mi) 



s + mi ' 
A^iiV2 |q||k| 



/s + mi 



-A. 



mi — s + ^k 



Sttv^ El + mi 

1 1,2 



J^4 ~ iVs — mi) 



■A2 + A3-A4 

\r_ i„|2 



fc2 



/s + mi 



-A, 



^7 



s — mi 
N1N2 |q| 
Sn^/s E2 + m2 

2fco 
qk 



Sny/s E2 + m2 
■^A. + A, 



Aa 



(mi - Ei)Ai 



fc2 



^/s - mi 



A, 



(2fco\/s — 3k ■ q) 
{2s-2ml-k^) ]A2 



^8 



2fc, 

+ (qoiVs-mi) - k ■ q)A3 

+ (k -q- qo{\/s- mi) + {Ei - mi){^/s + TO2))A4 
+ (gofc^ - kok ■ q)Az - {El - mi)(Vs + mi)Ae , 
^1^2 |k| 



(£;i +mi)Ai 



iny^ El + mi 



{2s-2mi-k^) ]A 



+ \^^{2koV^-3k-q) 

qk 

2ko 

+ {qaiVs + mi) - k ■ g) A3 

+ (k -q- qai^/s + mi) + {Ei + toi)(\/s - TO2))A4 
~ (gofc^ - kok ■ q)A5 - {El + mi)(Vs - mi)A6 , 



where 



N, = y^E^ + TO, , 



E^ 



1,2. 



We remark that, starting from eq. (jP.Sp . we have utilized 
the Lorentz condition k-e — 0. This is also valid when elec- 
troproduction is considered, since in that case the object 



is proportional to the electron-photon vertex of QED, 
and the condition then follows from current conservation. 

Restricting ourselves to s-and p-waves, we can use the 
CGLN-amplitudes to arrive at the multipoles Eq^, Afi+, 
etc.: 







Mi+ 




Mi- 


\-! 


\Ei+J 




and 






(Lo+ 


1 


Li+ 







dz 



( ^0 6 

4^1 ^4^2 -12^0,2 
^9-^1 9-^0 fi-^0,2 



1^0,2 \ 



12 



P0.2 1^A,3/ 



'Ti 
T2 




where Pi = Pi{z) are the usual Legendre polynomials and 
z is q • k, i.e. the cosine of the scattering angle in the cm. 
frame. Furthermore, the abbreviations 



Po,2 — Po ^ P2 



and 



were used. 

A word on units: Since we use 



1 



Pl,3 — Pi — P3 



Ana 



and normalize our Dirac spinors like uu = 2m, the in- 
variant amplitude -B5 has dimension GeV~^, as can be 
seen, 6.17., from the contribution of the graph correspond- 
ing to the 'Kroll-Ruderman' term. Therefore, the scatter- 
ing amplitude Tfi is dimensionless (see eq. (|D.1[) ). while 
the CGLN-amplitudes as well as the multipoles have di- 
mension GeV~^. 

The unpolarized differential cross section for —^ 
B M is given in terms of the CGLN-amplitudes as [50] 



1^2 P 



Re(J^iJ^4*) +Re(J^2,?^3*) 
(Re(,F3.F4*)-2Rc(,Fi.F2*)) cos tl 



;|,F4p 



.2' ' 2 
Re(JP3,^4*) 008^6* 



Re{TiTl + J-2^3) ) cos 



(D.5) 
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